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A new striking feature of hadron production in nuclear collisions is the large stopping of the 
participating nucleons in hadron-nucleus and nucleus-nucleus collisions. This enhanced baryon 
stopping can be understood introducing new diquark breaking mechanisms in multistring models of 
hadron production. Here we show, that similar diquark breaking mechanisms occur at high energy 
even in hadron-hadron collisions. This effect leads to significant changes in the extrapolation of 
, these models to Cosmic Ray energies. 
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\ I. INTRODUCTION 

<D 

A new feature of hadron production in nuclear collisions is the large stopping of the participating nucleons in 
' hadron-nucleus and nucleus-nucleus collisions. Experimental data demonstrating this effect have been presented in 

^ SB- 

Multistring fragmentation models like the Dual Parton Model (DPM) [|3|-|5[ or similar models did originally not 
show this enhanced stopping in nuclear collisions. Therefore, in order to incorporate the effect into multistring 
fragmentation models new diquark breaking DPM-diagrams acting in hadron-nucleus and nucleus-nucleus collisions 
were proposed by Capella and Kopeliovich fl an d investigated in detail by Capella and Collaborators Similar 
ideas were discussed by Vance and Gyulassy |10| and by Casado p3j] . 
^vq ' In the present paper we desribe a new diquark breaking diagram, which acts in nuclear collisions as well as in 
hadron-hadron collisions at sufficiently high energy. The new diagram becomes important only at energies well above 

■ the energies of the CERN-SPS heavy ion experiments which have been discussed in the papers quoted above. The 
effect of this new diagram modifies significantly the extrapolation of multistring fragmentation models to Cosmic Ray 

■ energies. This is our main interest in the present paper. 

The diquark breaking diagrams investigated by Capella and Collaborators |t],|8| use sea quarks, which in nuclear 
JL , collisions according to the Glauber model are needed to implement the multiple collisions in the nucleus. In the Dual 
Parton Model in the collider energy range there are further multiple collisions, again implemented using sea quarks at 
the ends of multiple chains, which occur due to the unitarization procedure. These are the sea quarks which provide 
!LJ ■ the diquark breaking in our new diagram. 
.£h ! To get definitive predictions from the old and the new diquark breaking diagrams, we introduce the diagrams 
into the Monte Carlo version of a multistring fragmentation model. DPMJET simulates hadron production in the 
framework of the Dual Parton Model with emphasis as well in the central as in the fragmentation region. 

Previous versions of the DPMJET event generator were described in detail Jl2|— [l5|| . Here we use the present version 
DPMJET-II.5 of the code. 

The most important applications of DPMJET so far were in Collaboration with Battistoni, Forti and others Jl^|2f|l 
for the simulation of the Cosmic Ray cascade in the HEMAS-DPM code system. 

Observations like rapidity plateaus and average transverse momenta rising with energy, KNO scaling violation, 
transverse momentum-multiplicity correlations and minijets pointed out, that soft and hard processes are closely re- 
lated. These properties were understood within the two-component Dual Parton Model f2^-|2g[|. The hard component 
is introduced applying lowest order of perturbative hard constituent scattering p9j . Single diffraction dissociation is 
represented by a triple-pomeron exchange (high mass single diffraction) and a low mass component. 

The Dual Parton model provides a framework not only for the study of hadron-hadron interactions, but also for 
the description of particle production in hadron-nucleus and nucleus-nucleus collisions at high energies. Within 
this model the high energy projectile undergoes a multiple scattering as formulated in Glaubers approach; particle 
production is realized by the fragmentation of colorless parton-parton chains constructed from the quark content of 
the interacting hadrons and nuclei. 

In Section II we describe the new diagrams introduced in multistring models like the DPM in order to get a better 
description of baryon stopping. In Section III we compare the model to data and in Section IV we discuss the 
properties of the model at the highest energies. 
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II. IMPLEMENTATION OF NEW DPM DIAGRAMS FOR AN IMPROVED DESCRIPTION OF 
BARYON STOPPING IN IHADRON-HADRON, HADRON-NUCLEUS AND NUCLEUS-NUCLEUS 

COLLISIONS 



A. Diquark fragmentation, the popcorn mechanism 



The fragmentation of diquarks is slightly more complicated than the fragmentation of a quark jet. As justified by 
Rossi and Veneziano |30[| the baryon can be pictured as made out of three quarks bound together by three strings 
which join in a so called string junction point. In diagrams we can characterize the baryons 

(i) by the three quarks and the string junction or 

(ii) by a quark and a diquark, in this case the string junction always goes with the diquark. 

In all the diagrams discussed in this section we will either plot the quark and the diquark or, if the diquark breaks, 
the three quarks and the string junction. The quarks and the two quarks of a diquark are plotted as solid lines, the 
string junction is plotted as a dashed line. 

There are two possibilities for the first fragmentation step of a diquark, see Fig]]]. Either we get in the first step 
a baryon, which contains both quarks of the diquark and the string junction or we get in the first step a meson 
containing only one of the two quarks and the baryon is produced in one of the following fragmentation step. This 
mechanism is well known, it is presented in the review on the Dual Parton Model 0] and it was investigated for in- 



stance in |3lJ,|32| . This mechanism was implemented from the beginning in the B AMJET-fragmentation code [|33 34 
used previously in DPMJET. This mechanism is also implemented under the name popcorn fragmentation in the Lund 
chain fragmentation model JETSET p5L|3 which is presently used in DPMJET. 




FIG. 1. (a) Conventional diquark fragmentation, the baryon is produced in the first fragmentation step, (b) The popcorn 
mechanism: a meson is produced in the first fragmentation step, the baryon appears in the second or later fragmentation steps. 
We plot the lines for the two quarks of the diquark as solid lines, the string junction is plotted as a dashed line. 

What happens in the model with the popcorn mechanism compared to the model without can be most easily 
seen looking at the proton rapidity distribution in p-p collisions. The two maxima in the target and projectile 
fragmentation region of the proton rapidity distribution shift by about half a unity to the center, these maxima 
become wider and correspondingly the dip in the center is reduced. At the same time the Feynman x distributions of 
mesons get a component at larger Feynman x. The effects in hadron-nuclcus and nucleus nucleus collisions arc quite 
similar. 

The popcorn effect as implemented in the JETSET fragmentation code |35| determines not only the fragmentation 
of a diquark as plottet in Fig. [l] A second effect of the popcorn mechanism is in baryon pair production inside 
the fragmenting chain. Here the popcorm mechanism determines the ratio of B — B ( baryon and antibaryon are 
neighbors in the chain) to B — M — B ( a meson is produced between the baryon and the antibaryon) production. 

The popcorn mechanism alone is not enough to explain the baryon stopping observed experimentally in hadron- 
nucleus and nucleus-nucleus collisions [0]|] , this will be discussed in more detail in Section III. 



B. New diquark breaking DPM diagrams in hadron— nucleus and nucleus— nucleus collisions 



New diquark breaking DPM-diagrams mainly of interest in hadron -nucleus and nucleus-nucleus collisions were 
proposed by Capella and Kopeliovich |6) and investigated in detail by Capella and Collaborators . Similar ideas 
were discussed by Vance and Gyulassy [Gil. Capella and Kopeliovich did discuss in detail their first diquark 
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breaking mechanism (see Fig.g, where this mechanism is characterized for nucleon-nucleon collisions), in this case the 
valence-diquark breaks into the two quarks, the baryon is produced in the second or in later fragmentation steps. 



c 



FIG. 2. The first C-K diquark breaking mechanism [g] plotted for a nucleon-nucleon collision. 

Looking only at the diagrams for a hadron-hadron collision, one does not sec any difference between the pop- 
corn mechanism and the first Capella-Kopeliovich (C-K) diquark breaking mechanism. However, the mechanism of 
this first (C-K) diquark breaking mechanism differs in detail and especially in hadron-nucleus and nucleus-nucleus 
collisions from the popcorn mechanism discussed above. 

Capella and Kopeliovich fj| write for the baryon rapidity distribution in the diquark breaking (DB) component 

^-(y)^[pM + p q A-y)} (2.1) 

where p Qv is the valence quark rapidity distribution. In hadron-nucleus (NA) collisions we have 



- <j dp + a DB [2.2) 



(J; 



where DP stands for diquark preserving. Studying the A behaviour of cr^g they find, that in N-A collisions cr^g 
rises faster with A than the diquark preserving component cr^p. The same happens in nucleus-nucleus collisions. 

Nevertheless, despite these important differences between the popcorn effect and the first C-K mechanism wc find 
that implementing the first C-K mechanism in DPMJET did not give any new feature of baryon stopping, which 
could not also be obtained from the popcorn mechanism. Therefore, we continue to use in DPMJET-II.5 the popcorn 
mechanism (which as discussed above also acts on baryon-antibaryon pair production in the middle of a chain) instead 
of the C-K mechanism, so far we did not find any argument to use both effects. 

More interesting is the second C-K mechanism, which was proposed but not discussed in detail in |J. This 
mechanism was discussed in detail by Capella and Collaborators [0-^] . In Fig.|] we plot first the diquark-conserving 
diagram for a nucleon-nucleus collisions with two participants of the target nucleus. This is the traditional way for 
such a collision in the DPM. 




FIG. 3. The diquark-conserving diagram for a nucleon-nucleus collision with two participants of the target nucleus. 

In Fig.[| we plot the second C-K diquark-breaking diagram for the same collision. Now the second valence quark 
from the broken diquark in the first C-K diagram in Fig.|] is replaced by a Glauber sea quark from the nucleon 
projectile. Therefore, we will call the mechanism the Glauber sea quark mechanism of baryon stopping GSQBS. The 
probability of such a diquark splitting rises if the considered nucleon is involved in more than two interactions. 
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Capella et al. JtJ] write for the rapidity distribution of the net baryon production AB = B — B in A-A collisions 

dN AA ^ AB , i n A ,dN£B _ dN AB 

— ^ — (y) = {y))n/n A + {n- n A ){ d y (y))n/n A \ ( 2 - 3 ) 

where n A is the average number of participants in each nucleus and n the average number of collisions. For n — h A 

dN &B 

only the diquark preserving component is present. — is the rapidity density of the diquark breaking component 

which behaves approximately like exp[— ^\vab — VMaxW- 

The GSQBS has been implemented in DPMJET-II.5 and we will see in Section III, that in nucleon-nucleus collisions 
and nucleus-nucleus collisions we are able with this mechanism to fill the dip in the baryon rapidity distributions at 
central rapidity in agreement to the experimental data. As discussed already in detail in this mechanism also 
contributes to increase the Hyperon production in nucleon-nucleus and nucleus-nucleus collisions. 
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FIG. 4. The Glauber sea quark mechanism of baryon stopping GSQBS for a nucleon-nucleus collision with two participants 
of the target nucleus. This is the second C-K diquark-breaking mechanism pi . 



C. The Casado diagram 

A new diagram was also introduced by Casado The diagram is plotted for a nucleon-nucleus collision with 
two participants of the target nucleus in Fig. |5[ Here the fragmenting diquark contains a Glauber sea quark and 
hadronizes like a valence diquark producing baryons mainly in the fragmentation region. However, the flavor content 
of the baryon is changed. This diagram gives another contribution to hyperon production in nuclear collisions. We 
have implemented the Casado diagram in DPMJET-II.5 and it is used in addition to the GSQBS diagram for all 
nuclear collisions compared in Section III to data. 




FIG. 5. The diagram introduced by Casado for a nucleon-nucleus collision with two participants of the target nucleus. 
The diquark contains one Glauber sea quark. 



D. The new diquark breaking diagram in hadron— hadron and nuclear collisions at high energy. A new 
extrapolations of DPMJET— II. 5 to Cosmic Ray energies. 

This new diagram follows necessarily from the GSQBS diagram if we go to high energy. At high energies we have 
multiple collisions even in hadron-hadron collisions due to the unitarization procedure. We call the sea quarks at the 
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ends of the additional chains in this case unitary sea quarks. The Glauber sea quarks are needed in nuclear collisions 
already at rather low energies, for instance at the energies of heavy ion collisions at the CERN-SPS. In contrast to 
this, unitary sea quarks appear in significant numbers in hadron-hadron and nuclear collisions only at rather high 
energies, for instance at the energies of the CERN-SPS collider or the TEVATRON collider, they are important in 
the Cosmic Ray energy region. 

With the unitary sea quarks at the ends of the chains from the secondary collisions we obtain a new mechanism for 
baryon stopping, which will become effective at very high energies. 




FIG. 6. Standard DPM diagram for a nucleon-nucleon interaction with one additional soft secondary interaction induced by 
the unitarzation procedure. There is one valence-valence and one sea-sea interaction, each represented by a pair of chains. 




FIG. 7. New DPM diagram for a nucleon-nucleon interaction with one additional soft secondary interaction induced by the 
unitarzation procedure. The diquark is split and an unitarity sea quark is used similar to Fig. ^ to shift the baryon in one of 
the chains. We call this the unitary sea quark mechanism for baryon stopping USQBS. 

In Fig.|| we plot the standard DPM diagram for a nucleon-nucleon interaction with two soft interactions induced 
by the unitarization procedure. There is one valence-valence and one sea-sea interaction, each represented by a pair 
of chains. In analogy to Fig. ^ we construct from this the new diagram for baryon stopping in Fig. |t]. The diquark 
is split and an unitary sea quark is used to have the baryon only in the second or later fragmentation steps in one of 
the chains. We call this the unitary sea quark mechanism for baryon stopping USQBS. Also here the probability for 
such a diquark splitting rises if there are more than 2 interactions of the hadrons involved. Obviously, this mechanism 
leads to Feynman x distributions of baryons in p-p collisions becoming softer and Feynman x distributions of mesons 
becoming harder than without the USQBS mechanism. 

One word about the implementation of the diquark breaking diagrams in DPMJET-II.5. In a Monte Carlo version 
of the Dual Parton Model we do of course not use the parametrizations used by Capella et al. In each event 

DPMJET first constructs the conventionel chains without diquark breaking. From this before implementing the 
chain fragmentation we construct the chain structure of the event including the diquark breaking mechanisms. In 
each elementary collision in DPMJET the number of participating Glauber sea quarks and unitary sea quarks is 
well known and given a probability for the diquark breaking diagram, each possible sea quark gets the chance to get 
involved in the diquark breaking mechanism. Also in DPMJET we do not need any parametrization for the rapidity 
distributions of the DP and DB components, with given x-values of the participating diquarks, valence quarks and 
sea quarks and using the JETSET chain fragmentation we get the resulting hadrons at the correct positions, energy 
conservation in the event is satisfied in this way. If the diquark is broken, we have to convert the x of the diquark 
into the two x— values of the two valence quarks. At low energy and for original diquark-quark chains of low invariant 
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mass we have to reject the diquark breaking at times. If the diquark breaking is kinematically impossible, we retain 
the original DP mechanism. 

In contrast to the GSQBS mechanism, which leads already to effects in nuclear collisions at the energy of the CERN- 
SPS, we have at present no data available to prove that this USQBS mechanism is a valid extension of the DPM. There 
is no data for baryon stopping in proton-proton or antiproton-proton collisions at collider energies. Unfortunately, 
the fragmentation region at large Feynman x has not been investigated experimentally with enough detail at any of 
the hadron-hadron colliders. However, if the GSQBS mechanism is the correct mechanism responsable for the baryon 
stopping effects found in nuclear collisions, then also the USQBS mechanism should modify the collisions at collider 
energies and beyond. In Section IV, where we present the properties of DPMJET-II.5 at Cosmic Ray energies we 
present always two extrapolations of the model. In one version (shortly characterized as version 55 in the plots) we 
have both the GSQBS and USQBS mechanism, in a second version (characterized in the plots as version 50) we have 
only the GSQBS mechanism. Only collider experiments on baryon stopping or Cosmic Ray observations can prove, 
which version is the better one. From theoretical prejudices however there would be the claim, that the version with 
both the GSQBS and USQBS mechanisms is to be prefered. 

The version of the model with both GSQBS and USQBS (designated by 55 in the plots) is the new version of the 
model with all new baryon stopping effects. The model with only GSQBS (designated by 50 in the plots) serves as 
reference only, it corresponds in its high energy behaviour closely to previous versions of DPMJET or to other models 
without the new baryon stopping mechanisms. 

E. New parameters connected with the diquark breaking diagrams 

For each of the new diquark breaking diagrams described in this Section we have to introduce a new parameter. 
These parameters give the probability for the diquark breaking mechanisms to occur, given a suitable sea quark is 
available and given that the diquark breaking mechanism is kinematically allowed. For an original diquark-quark 
chain of small invariant mass, which originally just fragments into two hadrons, the diquark breaking is often not 
allowed at small energies. 

The optimum values of the new parameters are determined by comparing DPMJET-II.5 with experimental data 
on net-baryon distibutions. We use for the GSQBS and USQBS mechanisms the default parameters 0.45 and for the 
Casado diagram the default parameter 0.5. 

F. Production of strange particles 

Capella et al 0-^) discuss the new diagrams for baryon stopping also with the interest in strangeness production 
in nuclear collision. They discuss in detail the production of hyperons and antihyperons. We will not repeat their 
discussion here, we only shortly summarize the different mechanisms, which contribute to strangeness production. 

Studies of strangeness production within DPMJET were given in |37] f39|] . Enhanced generation of strange particles, 
in particular of strange antibaryons, has been proposed as a signal for the formation of quark gluon plasma in dense 
hadronic matter [ ^0|j4l] ]. Recent data from experiments at the CERN SPS have already been interpreted within this 
scheme. However, we find it worthwhile to pursue the study of tested conventional models without QGP formation 
like the DPM before drawing final conclusions. The DPM is an independent string model. Since the individual 
strings are universal building blocks of the model, the ratio of produced strange particles over non-strange ones will 
be approximately the same in all reactions. However, since some strings contain sea quarks at one or both ends and 
since strange quarks are present in the proton sea, it is clear, that, by increasing the number of those strings, the ratio 
of strange over non-strange particles will increase. This will be the case for instance, when increasing the centrality 
in a nucleus nucleus collision. It is obvious, that the numerical importance of the effect will depend on the assumed 
fraction of strange over non-strange quarks in the proton sea. The rather extreme case leading to a maximum increase 
of strangeness would be to assume a SU(3) symmetric sea (equal numbers of it, d and s flavors). We express the 
amount of SU(3) symmetry of the sea chain ends by our parameter s sea defined as s sea = 2 < s s >/(< u s > + < d s >) 
where the < q s > give the average numbers of sea quarks at the sea chain ends. Usually, DPMJET uses the default 
value s sea = 0.5. 

There are a number of effects, which change the number of strange hadrons, especially strange and multistrange 
baryons and antibaryons in the model: 

(i) The presence of sea qq and qq diquarks at sea-chain ends p2LE3| . 

(ii) The popcorn effect. 

(iii) The diquark breaking diagrams discussed in Section II. C, see also [H-pl. 
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(iv) The secondary interactions of co-moving produced particles H,f45[], see also Section II. G. 

(v) Effects of string fusion and percolation of strings |48|j4g| 1 . 

It is certainly beyond the scope of the present paper, to discuss all of these effects in detail. The effects (i) to (iv) 
are contained in DPMJET-II.5. 

G. Final state interactions of co— moving secondaries in nuclear collisions 

Analyzing the rapidity distributions of produced A and A perticles in central heavy ion collisions at CERN-SPS 
energies within the Dual Parton Model Capella et al Jl4|,fl5| noted the need for secondary interactions of co-moving 
secondaries to understand the data. Secondary interactions of produced hadrons are also considered in other Monte 
Carlo models for hadron production in nuclear collisions [pp|-p2[ . 

Capella et al @,f|5| only introduced the following secondary ir + N — * K + A and 7r + iV — > K + A. A reasonable 
cross section for these reactions of a » 1.5 mb was needed in 0^^] to understand the data. All of this was done in 
|Il| . [l5f using analytical methods, not a Monte Carlo event generator. Final state interactions of co-moving secondaries 
was also used in |5^,|54|] to explain the J/ip supression in Pb-Pb collisions. 

The method of Capella et al @,|^| was implemented in DPMJET already in 1995. The only reason to discuss 
final state interactions of co-moving secondaries here is our comparison in the following Section of DPMJET-II.5 in 
this version with experimental data. In DPMJET the method is used to modify the Monte Carlo events and we use 
only a cross section of a w 1.0 mb. First results reproducing essentially the results of Jl4|,[l5| are described shortly 
in [T|] and in more detail in a unpublished code write-up p4[ . We would like to stress, this method as implemented 
with only the two reactions given above, can only be considered as a first preliminary step. A better method of final 
state interactions should be implemented as a Monte Carlo method from the beginning and more types of secondary 
interactions should be taken into account including the reverse reactions. With the present method we can not expect 
to obtain reasonable results in situations with much higher secondary particle densities than in the heavy ion collisions 
at the CERN-SPS energy. Therefore, it is not recommended to use the secondary interaction option of DPMJET at 
RHIC or CERN-LHC energies. 

III. COMPARING DPMJET-II.5 TO DATA WITH EMPHASIS TO BARYON STOPPING 
A. Comparing to data on leading particle production in hadron— hadron collisions 

The leading particle production is very important for the Cosmic Ray cascade simulation. Also, the leading baryons 
are directly influenced by the diquark breaking mechanisms. 

In Fig. |8| we compare the model with NA-22 data on the Feynman-x distribution of positively charged hadrons 
produced in 250 GeV pp collisions. In Fig. || we observe in particular (compared to previous DPMJET versions) a 
better agreement of the model predictions with the forward production of protons. 

Distributions of leading protons are compared more directly to data in Fig. ^. The Figure is adopted from talks of 



Engel 1 55 59]. But of course, we present now the DPMJET-II.5 results. In Fig. O we compare the distribution in the 



energy fraction xi a b carried by the leading proton. The data are photoproduction and DIS measurements from the 



HERA Collider at «/s = 200 GeV J57J5J. We compare to DPMJET-II.5 for p-p collisions at ^ = 200 GeV. The 
forward production of leading protons is not expected to depend strongly on the reaction channel. It is found, that 
DPMJET-II.5 agrees much better to the data than older DPMJET versions, see J5(|. This comparison demonstrates, 
that the leading proton distribution at y/s = 200 GeV is still rather flat. As shown in Fig. || at ^/s = 200 GeV 
both versions of DPMJET-II.5 (discussed in Section II) with GSQBS and USQBS (designated in the plot by 55) and 
with only GSQBS (designated in the plot by 50), lead still to very similar xi a b distributions of secondary protons. 
Unfortunately therefore, these data cannot be used to discriminate between the two DPMJET-II.5 versions. 

In Fig. n(| we compare the leading baryon distribution in the energy fraction xi a b according to both versions of 
DPMJET-II.5 at = 200 GeV and 100 TeV. At */s = 200 GeV we find at x lab obove 0.5 like for the leading proton 
distribution in Fig. ^ hardly any difference between both versions of the model. At s/s = 100 TeV the differences 
between both versions are significant even for xi a b above 0.5. At smaller xi a b values, say below x^ =0.2 we find 
significant differences between both versions at both energies. The xi a b distribution of leading protons looks rather 
similar to Fig. [l^. We conclude leading proton or leading baryon distributions in the energy fraction xi a b below xi a b 
=0.5 could be quite usefull to demonstrate the effect of the USQBS diagrams in experimental data. 
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FIG. 8. Comparison of Feynman-x distributions of positively char ged hadrons produced in proton-proton collisions at 250 
GeV. The experimental data are from the EHS-NA22 Collaboration [p9| . 
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FIG. 9. Energy fraction xi a b carried by the leading proton. The data are photoproduction and DIS measurements at y/s = 
200 GeV [^7|j5^] shown as symbols compared to both versions of DPMJET-II.5 for p-p collisions at y/s = 200 GeV. 
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FIG. 10. Energy fraction xi a b carried by the leading baryon. at y/s =200 GeV and 100 TeV according to both versions of 
DPMJET-II.5 for p-p collisions. 



B. Comparing to data in hadron— nucleus and nucleus— nucleus collisions 



We now turn to collisions with nuclei. We want to compare the model with the diquark breaking mechanisms to 
data in hadron-nucleus and nucleus-nucleus collisions. Therefore we have to demonstrate thet the model reproduces 
correctly the global properties of these nuclear collisions. More comparisons are presented in |l6] |. 

In Table 1 we compare average multiplicities of negatively charged hadrons calculated with DPMJET-II.5 in 
hadron-hadron, minimum bias hadron-nucleus and central nucleus-nucleus collisions with experimental data. 

Table 1. Comparison of average multiplicities of produced negatively charged hadrons in proton-proton, minimum 
bias proton-nucleus and central nucleus-nucleus collisions at 200 GeV. 



Collision 


DPMJET-II.5 


Exp. 


Reference 


p-p 


2.85 


2.85 ± 0.03 


go 


1 


p-S 


5.10 


5.0 ±0.2 


H 




p-Ar 


5.30 


5.39 ±0.17 


n 




P-Ag 


6.18 


6.2 ±0.2 


d 




p-Xe 


6.43 


6.84 ±0.13 


1 




p-Au 


6.81 


7.0 ±0.4 


n 




p-Au 


6.81 


7.3 ±0.3 






S-S central 


103 


98 ±3 


1 


1 




S-Ag central 


174 


186 ±11 


1 


1 




S-Au central 


202 


225 ±12 


I 


1 









Table 2. Comparison of average multiplicities of produced strange hadrons in central S-S and S-Ag collisions 
at 200 GeV. For DPMJET-II.5 we give the results for the model with secondary interactions of comovers and it is 
assumed that E° and S° have decayed. The experimental data are from the NA35 Collaboration p5| . 



Reaction 


Particle 


DPMJET-II.5 with sec. int. 


Exp. Q 


S-S 


A 


7.3 


9.4 ± 1.0 


S-S 


A 


1.22 


2.6 ± 0.3 


S-S 


s 


9.1 


10.5 ± 1.7 


S-S 


K+ 


11.4 


12.5 ± 0.4 


S-S 


K- 


6.9 


6.9 ± 0.4 


S-Ag 


A 


13.1 


15.2 ±1.2 


S-Ag 


A 


2.0 


2.6 ± 0.3 


S-Ag 




15.5 


15.5 ±1.5 


S-Ag 


K+ 


19.7 


17.4 ± 1.0 


S-Ag 


K~ 


11.4 


9.6 ± 1.0 



In Table 2 we give the average multiplicities of produced strange hadrons in central S-S and S-Ag collisions at 200 
GeV. We give the DPMJET-II.5 results for the model with secondary interactions and it is assumed, that S° and S° 
have decayed. The experimental data for central S-S and S-Ag collisions are from the NA35 Collaboration [ jS5| . 

C. Baryon stopping 

We present first the DPMJET predictions at pi a i, = 200 GeV/c for net baryon rapidity distributions in the original 
model without the new GSQBS and USQBS diagrams. In Fig. [ll] we present the leading net proton (p — p) rapidity 
distribution dN p / dy — dNp j dy in p-p, p-S and central S-S collisions. In Fig. [l2]we present the net A (A — A) rapidity 
distribution dN^/dy — dN^/dy in p-p, p-S and central S-S collisions. In p-p collisions, which at the given collision 
energy are hardly modified by the new USQBS diagram, we observe a dip at central rapidity. This dip is also present 
in DPMJET without the new GSQBS and USQBS diagrams in p-S and central A-A collisions. This disagrees to the 
data presented in the next Figures. 



10 



10 



- A 



1 



i 1 r 

p+p u-A-J ■ 

P+S r-A-J 

S+S central ;■•-+- 



++ + +++ +++++ , , +++++++ ++ +++ + 



+ A 

A 

A A A 

0.1 fcr 



+■ 



aAaA4S a 



J_ 



J I L 







1 



5 



2 3 4 
Rapidity j/ 

FIG. 11. Net proton (p — p) rapidity distribution in p-p, p-S and central S-S. collisions. Calculated with DPMJET-II.5 
without the new diagrams modifying baryon stopping. 
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FIG. 12. Net A (A — A) rapidity distribution in p-p, p-S and central S-S. collisions. Calculated with DPMJET-II.5 without 
the new GSQBS and USQBS diagrams modifying baryon stopping. 

In Fig. (l| and [jlj we compare the net-proton distributions according to the full DPMJET II. 5 model with data in 
p-S and p-Au collisions Now the dips at central rapidity are filled in the model, we observe like in the data at 
central rapidity a flat net-proton rapidity distribution. 
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FIG. 13. Net proton (p — p) rapidity distribution in p-S collisions. The DPMJET-II.5 results are compared with data 
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FIG. 14. Net proton (p — p) rapidity distribution in p-Au collisions. The DPMJET-II.5 results are compared with data 



In Fig. [15| we compare the full DPMJET II. 5 model in the version with secondary interactions with data on 
net-proton production in central S-S collisions. We give the DPMJET results with and without spectator protons 
(evaporation protons from the residual nuclei). Also here the dip at central rapidity in the model has disappeared 
(compare to Fig.|ll|), however, the agreement to the data S is not perfect. There is a significant disagreement in the 
fragmentation regions of the two nuclei, especially with the DPMJET version including the spectator protons. The 
spectator evaporation protons are apparently not included in the data. 
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FIG. 15. Net proton (p — p) rapidity distribution in central S-S collisions. The DPMJET-II.5 results are compared with 
data Q. DPMJET-II.5 is used with secondary interactions of comovers and we present the DPMJET results with and without 
proton spectators. 



In Fig. |T| we compare the full model to data on net-A production in p-S collisions. In FigJlq the dip at central 
rapidity has disappeared in the model. 
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FIG. 16. Net A (A — A) rapidity distribution in p-S collisions. The DPMJET-II.5 results are compared with data 0|. 

IV. PROPERTIES OF THE MODEL IN THE HIGHEST ENERGY REGION, THE CONSEQUENCES OF 
THE BARYON STOPPING MECHANISMS AT COSMIC RAY ENERGIES. 



In all plots in this Section we present the DPMJET-II.5 results for the two versions of the model discussed already 
in Section II. C: 

(i) With only the new diagram for baryon stopping GSQBS, this version is characterized in the plots by 50, this 
version corresponds in its high energy behaviour approximately to the previous version of the model DPMJET-II.4. 

(ii) With both new diagrams for baryon stopping GSQBS and USQBS, this version is characterized in the plots by 
55, the hight energy behaviour of this version is new, at high energies the baryons carry 5 to 10 percent less energy 
than in version 50, correspondingly the mesons carry more energy than in version 50. 

We start as example with one plot where the difference between the two versions of the model is rather insignificant. 
With rising energy the fraction of strange hadrons rises slowly. In Fig. [l7| t he K/7T ratio according to DPMJET-II.5 
for pp or pp collisions is compared to data from the E735 Collaboration pgfl. 
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We observe similar insignificant differences between both versions of the model in plots of the average transverse 
momentum < p± > c h as function of the collision energy or of the average charged multiplicity < n c h > as function of 
the collision energy. 
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FIG. 17. K/-7T ratios in pp or pp collisions as function of the cms energy t/s. The DPMJET-II.5 calculation is compared with 
data collected from the E735 Collaboration at Fermilab |6^ |. 

Next we discuss plots with significant differences between the versions with (55 in the plots) and without (50 in the 
plots) the USQBS mechanism. 

Following for instance the basic discussion of ||S7| , we introduce a variable xi a b similarly to Feynman-XF , but this 
time in the lab-frame : 



Xlab 



E 



(4.1) 



Ei is the lab-energy of a secondary particle i and E is the lab-energy of the projectile in a h-nucleus collision. 
We introduce xi a b distributions F{xi a b) ■ 
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FIG. 18. Spectrum weighted moments for pion production in p-N collisions as function of the (nucleon-nucleon) cms energy 



Fi(Xlab) = Xlab 



dxi ab 



(4.2) 
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We note that the Feynman-x distribution at positive in the projectile fragmentation region is a very good 
approximation to the xi a b distribution. 

The cosmic ray spectrum-weighted moments in p-A collisions are now defined as moments of the F(xi a b) '■ 

ff- A = i\xi ab y- 1 Ff- A {x lab )dx lab (4.3) 
Jo 

Here —7 ~ -1.7 is the power of the integral cosmic ray energy spectrum and A represents the target nucleus. 



0.013 
1.0125 



+" 



DPMJET-25 p-N 55 -•- 
DPMJET-25 p-N 50 




lL 



J U I I L 



0.1 1 10 100 1000 

Energy <Js [TeV] 

FIG. 19. Spectrum weighted moments for Kaon production in p-N collisions as function of the (nucleon-nucleon) cms energy 



The spectrum- weighted moments for nucleon-air collisions, as discussed in 
of energetic particles in the atmosphere. 

We also introduce the energy fraction Kf~ : 



determine the uncorrelated fluxes 



K. 



p-A 



F[ A (xiab)dXl ab 



(4.4) 



As for xi a b, the upper limit for K is 1 in h nucleus collisions. 

In FigsjlS] we present the spectrum weighted moments summed over pions of both charges in p-N collisions as 
function of the cms energy ^fs per nucleon. In Figsjl^ the moments are given for charged Kaon production also in 
p-N collisions. 

In Fig.^0] and [2l] we present again for pp anf p-N collisions the energy fractions K for net baryons B — B (baryon 
minus antibaryon) and charged pion production. The energy fraction K B _g is always smaller than the energy fraction 
of baryons Kb- The difference between both is the energy fraction going into antibaryons Kg which is equal to the 
energy fraction carried by the baryons which are produced in baryon-antibaryon pairs. 
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FIG. 20. Laboratory energy fractions for net baryons (baryon minus antibaryon) B — B, baryons B and charged pion 
production in p-p collisions as function of the (nucleon-nucleon) cms energy ^/s. 
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FIG. 21. Laboratory energy fractions for net baryons (baryon minus antibaryon) B — B, baryons B and charged pion 
production in p-N collisions as function of the (nucleon-nucleon) cms energy -y/s. 



We find in DPMJET-II.5 all average values characterizing hadron production: the average transverse momenta, 



the charged multiplicities and the moments in Figs. 18, [19| [20| and |2l| to change smoothly with energy in most cases 
just like the logarithm of the energy. 

So far, we have not found any experimental data to favour either the versions of DPMJET-II.5 (i or 50) with only 
GSQBS or (ii or 55) with GSQBS and USQBS. Clearly, the version prefered by theoretical prejudices is (ii or 55). In 
this version we have a better Feynman scaling of meson distibutions and spectrum weighted moments and we have a 
faster decrease with the collision energy of the energy fractions into secondary baryons. 



V. SUMMARY 



In the present paper we discuss new diquark breaking mechanisms, which lead to a better agreement with experi- 
mental data on baryon stopping in hadron-nucleus and nucleus-nucleus collisions in multichain fragmentation models. 
This is demonstrated using a particular Monte Carlo version DPMJET-II.5 of the Dual Parton Model. A new diquark 



1G 



breaking diagram USQBS is emphasized, which is unimportant at the energies of the present baryon stopping data. 
This new diquark breaking diagram becomes very important at high energies and it changes significantly the extrapo- 
lation of the multistring models to Cosmic Ray energies. At the same time the new diquark breaking mechanisms do 
not spoil the good agreement of the model to leading proton Feynman x distributions at present accelerator energies. 
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